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CsJ ■ Abstract 

The production of rj' mesons in the reactions pp — > ppr]' and pn — > pnij' at threshold-near 

energies is analyzed within a covariant effective meson-nucleon theory. The description of cross 

O ' section and angular distributions of the available data in this kinematical region in the pp channel 

G^ ' is accomplished by including meson currents and nucleon currents with the resonances 5ii(1650), 

Pii(1710) and Pi3(1720). Predictions for the pn channel are given. The di-electron production 

Q^ I from subsequent rj' Dalitz decay r]' -^ 77* -^ je'^e is also calculated and numerical results are 

^H ! 

^D ■ presented for intermediate energy and kinematics of possible experiments with HADES, CLAS 
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t:^' : and KEK-PS. 
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I. INTRODUCTION 

The pseudo-scalar mesons 77 and rj represent a subject of considerable interest since 
some time which has been addressed in many investigations (see, e.g. [1-5]). The physical 
states of ?7, r( mesons can be constructed as a superposition of r^o, f]^ members of a SU(3) 
pseudoscalar nonet with one mixing angle. If the coupling constants (defined usually 
by the relation (0|J"5|P(p)) = ifpPfi, where P is a member of the SU(3) nonet, and 
a = — 8) follow the same pattern as the state mixing then they must be related to 
each other as [6] /^ = /8Cos6', /° = — /osin6', /^, = /gsin^, /°, = /qcos^, where 9 
is the state mixing angle. A combined phenomenological analysis [7] of data on two- 
photon decays of rj and rj' as well as 777 transition form factors shows indeed that one can 
parameterize the relation among the coupling constants by such a relation, however, with 
an angle essentially differing from the state mixing angle. The interrelation between the 
two angles and different decay constants can be found by considering the divergences of 
axial vector currents including the U{1) axial anomaly, which couples the f]' state with 
gluon fields [8, 9]. 

The U{1) axial anomaly plays also a crucial role in understanding the physical masses 
of the 1] and rj' mesons from the point of view of spontaneous breaking of chiral symmetry. 
It is known [3, 4] that the non-zero divergence of the axial current J^g in the chiral limit 
in presence of interactions is due to couplings with gluons arising after renormalization. 
This anomaly is responsible for the generation of the heavy mass of the physical f]' meson. 
Yet, the anomaly is connected with a violation of the OZI rule for the nucleon and seems 
also to point to a small strangeness admixture in the proton wave function [10]. 

This also has a direct consequence in deep-inelastic scattering. For instance, the 
Goldberger-Treiman relation (see, e.g., [10] and references therein) relates the nucleon- 
nucleon-r]' coupling constant gNNrf (which is non-zero because of OZI rule violation) with 
the fiavor-singlet axial constant G'a(O) and to the coupling g^Ng with gluons. If the gluon 
contribution is ignored in this relation then the EMC measurements lead to the so-called 
"spin crisis" [11]. A knowledge of all terms in the Goldberger-Treiman would allow to 
understand the origin of the spin crisis. However, neither g^Nr]' nor g^Ng have been mea- 



sured directly in experiments. From this point of view an investigation of the pp -^ pprj' 
reactions can be considered a tool for supplying additional information on a wide range 
of fundamental theoretical issues. 

The strong decay channels of 7] and rj' are highly suppressed by various constraints 
(such as G parity) so that their widths are extremely narrow implying a relatively long 
life-time and, consequently, a suitable means to investigate different theoretical aspects 
of the OZI rule, U{1) axial anomaly, role of gluons in the NNt]' vertex etc. A further 
remarkable fact is that near the threshold the invariant mass of the NNt]' system in 
corresponding reactions is in the region of heavy nucleon resonances, i.e. resonances with 
isospin i can be investigated via these processes. Furthermore, resonances weakly bound 
with nucleons (the so-called "missing resonances" [12, 13]) can be studied in principle. 
First SAPHIR data on rj' photo-production [14] has shown a strong rise of the cross section 
near the threshold which indicates a possibly large contribution from nucleon resonance 
excitation. Moreover, if one assumes that the production mechanism is solely governed by 
resonances then the new states 5*11(1897) and Pii(1986) can be predicted [14]. However, 
further analysis of the SAPHIR and CLAS [15, 16] data has shown [17-19] that data are 
also compatible with a larger number of resonances, including higher spin states | and |, 
i.e., data are still too scarce to firmly determine the properties of the relevant resonances. 

Another aspect of rj' production in elementary hadron reactions is that the subsequent 
Dalitz decay may constitute a prominent source of di-electrons in intermediate-energy 
heavy-ion collisions. Indeed, the recent HADES data [20] exhibit a sizeable di-electron 
yield at large invariant masses which, besides the contribution from the vector mesons p 
and uj [21], can be due to rj/r]' mesons too. The rj Dalitz decay has been quantified [21] 
as prominent contribution to the invariant mass spectrum, indeed. One of the primary 
aims of the HADES experiments [22] is to seek for signal of chiral symmetry restoration 
in compressed nuclear matter. For such an endeavor one needs a good control of the 
background processes, including the f]' Dalitz decay, in particular at higher beam energies, 
as becoming accessible at SISlOO within the FAIR project [23]. 

The f]' Dalitz decays depend on the pseudo-scalar transition form factor, which encodes 
hadronic information accessible in first-principle QCD calculations or QCD sum rules [5]. 



The Dalitz decay process of a pseudo-scalar meson P can be presented as P — ;► 7 + 7* — ;► 
7 + e^ + e+. Obviously, the probability of emitting a virtual photon is governed by the 
dynamical electromagnetic structure of the "dressed" transition vertex P ^ 7* which 
is encoded in the transition form factors [5, 24-27]. If the decaying particle were point 
like, then calculations of 7* mass distributions and decay widths would be straightforward 
along standard quantum electrodynamics (QED) techniques. Deviations of the measured 
quantities from the QED predictions directly reflect the effects of the form factors and, 
thus the internal hadron structure, and, consequently, can serve as experimental tests to 
discriminate the different theoretical models in the non-perturbative QCD regime. Yet, 
calculations of the transition form factor within perturbative QCD involve triangle and 
box diagrams, which are tightly connected to the f/(l) axial anomaly so that form factors 
provide additional information on the nature of anomalies in quantum field theories, here 
QCD. 

In the present paper we study the di-electron production from Dalitz decay of 77' mesons 
in pp and pn reactions at beam energies of a few GeV for kinematical conditions corre- 
sponding to the HADES setup [22]. Our focus is to elaborate a model which provides 
a reliable energy dependence of the total cross section of rj' meson production and its 
subsequent decay into a di-electron and a photon at threshold-near beam energies. To 
this end we calculate the dependence of the differential cross section for the reaction 
pp — > pp'~fe~^e~ upon the invariant mass of the subsystem 7e"'"e~ around the pole masses 
of 1]'. We calculate the invariant mass distribution of di-electrons in a suitable kinematical 
range as a function of the di-electron invariant mass and argue that such a quantity, nor- 
malized to the real photon point and supplemented by some specific kinematical factor, 
represents the transition form factor. In such way the extraction of the transition form 
factor becomes accessible. 

1]' production processes have been analyzed in several papers. In Refs. [17-19, 28], 
photo-production of rj' has been studied to elucidate the role of resonances. An investi- 
gation of the 1]' production in A^A^ reactions has been performed in Refs. [29, 30], and a 
consistent combined analysis of rj' production in the reactions pp -^ pprj' and 7p -^ prj' 
has been attempted in Ref. [31]. 



Our paper is organized as follows. Section II is devoted to the theoretical background 
for dealing with the reaction pp -^ pprj' -^ pp^e^e~ . It is essentially based on the effective 
model [26, 27, 32, 33] proposed to describe bremsstrahlung and vector meson production 
in A^A^ reactions. The model is based on direct calculations of the relevant Feynman dia- 
grams within a covariant phenomenological meson-nucleon theory. The model parameters 
have been fixed from independent experiments and adjusted to achieve a good descrip- 
tion of the available experimental data. In the present paper, diagrams with excitation 
of nucleon resonances with masses close to the masses of a nucleon plus rj' meson have 
been also included. These are 5*11(1650), Pii(1710) and Pi3(1720) resonances. The cor- 
responding effective constants, whenever possible, have been obtained from the known 
decay widths of direct decay into the rj' channel or radiative decay with subsequent use of 
vector meson dominance conjecture. Also we use effective constants commonly adopted 
in the literature and obtained from different considerations, e.g., SU(3) symmetry, fit of 
photo-absorbtion reaction etc. [29]. In Section III, total cross sections and angular dis- 
tributions for 7]' production in A^A^ reactions are presented. Comparison with available 
experimental data is also performed. In Section IV, the rj' Dalitz decay is considered. 
Results of calculations of the invariant mass distribution as a function of the di-electron 
mass are presented for different initial energies. The role of the transition form factor in 
such processes is investigated as well. Our conclusions are summarized in Section V. 

II. THE MODEL 

We consider the reaction 

Ni + N2-^Ni + N2 + 'n' ^Ni + N2 + i + e^ + e-. (2.1) 

The invariant cross section is 

where the factor l/n\ accounts for n identical particles in the final state, |Tp denotes the 
invariant amplitude squared, m is the nucleon mass, s denotes the invariant mass of a 



particle or a sub-system of particles with total four-momentum p, i.e., s = p^, and drf 
is the invariant phase volume. The kinematical factor A, describing the incident flux, is 
defined as A(a:^, y'^, z^) = {x^ — {y + zy){x^ — iy — z)"^)- Eventually, P[, P2, h, k2, k^ are 
the momenta of the final nucleons, leptons and photon, respectively. 

A. Interaction Lagrangians 

The invariant amplitude T is evaluated here within a phenomenological meson-nucleon 
theory based on effective interaction Lagrangians (see [26, 27, 29, 32, 33]) which include 
(i) isoscalar mesons: scalar (a), pseudo-scalar {1]), vector (u;), (ii) isovector mesons: scalar 
(oo), pseudo-scalar (tt) and vector (p): 

C^NN = gaNNNN^^, (2.3) 

J^aoNN = gaoNNN{T^ao)N, (2.4) 

C^NN = -^-^^Ni,Yd,{T^^)N, (2.5) 

C,NN = -^^^^N-f,Ydf.%N, (2.6) 

CpNN = -9pNN [ni.tN^p^ - l^iVa^.riV^'^*/) , (2.7) 

C^NN = -g.NN {N-fpN^i:; - ^Nap.Nd'^^!:^^ , (2.8) 

where N and ^ denote the nucleon and meson fields, respectively, and bold face letters 
stand for iso-vectors. All couplings with off-mass shell particles are dressed by monopole 
form factors Fm = (A|,j — fij^) / (A|^ — A;|^), where kj^ is the 4-momentum of a virtual 
particle with mass fiM- 

To describe the Dalitz decay of the r]' meson Eqs. (2.3)-(2.8) must be supplemented 
with the corresponding Lagrangians describing the interaction of the electromagnetic field 
Ap with di-electrons (// = e+e^) and with the rj' meson, 

C.,ii = ~e{^ap^i)A>', (2.9) 

C^,^^ = f^,^^[eapp.d^A''d''A^)<^^,, (2.10) 

where £0123 = —1- Lagrangians (2.3)-(2.10) determine the S matrix from which one 



generates the corresponding tree-level Feynman diagrams describing the " bremsstrahlung" 
of 1]' off nucleons and the subsequent Dalitz decay rj' -^ 7e+e^. Usually such diagrams 
are called in the literature as nucleonic currents in Dalitz decay of mesons. 

The f]' meson can be produced also by an internal conversion of the exchanged mesons, 
the so-called conversion or meson current. The dominant exchange mesons in this case 
are u and p mesons with the interaction Lagrangians 



9r]'pp 



^V- = -f^^M-/3 (^"C5"€) '^v', (2-11) 



^v'PP = -|f ^M-/3 [d'^;d^^'p) K- (2-12) 

The corresponding diagrams are the same as in rj production [27], see Fig. 1. 

In the threshold- near kinematics for the rj' production in NN reactions there are several 
firmly determined nucleon resonances (four stars, according to the Particle Data Group 
classification [34]) and a large number of the so-called "missing resonances" [12, 13], which, 
in principle, can be excited in these processes. Consequently, rj' production can serve as 
additional tool to search for and to investigate the missing resonances in the mass interval 
1.5 — 2.0 GeV/c^. As mentioned in the Introduction, in Refs. [18, 28, 30] a preliminary 
analysis of SAPHIR data [14] has been performed assuming that only resonance currents 
contribute in this kinematical region. Then, excitations of missing 6*11 and Pu resonances 
have been considered and their pole masses and widths are estimated. A further combined 
analysis [19], including also the CLAS data [15, 16], has demonstrated that excitations 
of only two resonances are not sufficient to describe the data, and at least two more 
missing resonances, P13 and D13, ought to be involved into the calculations. Eventually, 
a systematic analysis of the rj' production in photo and A^A^ reactions [31] has shown 
that an equally well description of data could be achieved with several diverse sets of 
diagrams, which include different numbers of known and missing resonances. This means 
that up to now the available data is too scarce to determine unambiguously the kind 
and characteristics of resonances contributing in this kinematical region. Moreover, since 
in photo and A^A^ reactions one can observe excitations of different resonances, it is not 
mandatory to analyze simultaneously photo and A^A^ data within the same set of diagrams. 
We are interested in finding reliable parameterizations of the energy dependence of the 



total cross section and angular distributions near the threshold. Consequently, in order to 
reduce the number of free parameters and to avoid additional ambiguities, in the present 
analysis we consider only the known (four stars) lowest spin resonances 5*11(1650) with odd 
parity, and Pii(1710) and Pi3(1720) with even parity. The corresponding nucleon-meson- 
resonance interaction Lagrangians can be found in Ref. [29]. It should be mentioned 
that the inclusion of higher spin resonances leads to additional uncertainties. It is known 
that the Lagrangian for particles with spins s > | possesses an additional symmetry, 
i.e., additional free parameters (see Ref. [35] and references therein quoted). Thus, the 
invariance of the Lagrangian with respect to the point transformation causes an additional 
transformation parameter A and the of-mass shell parameter z. In the present calculation 
we adopt A = —1 and z = —1/2. Also, the choice of the form of higher-spin propagators 
has been a subject of discussion in the literature [33, 36-39] with respect to the choice of 
the spin projector operator Ps (off mass shell or on mass shell) and the order to write the 
product of the energy projection operator P/v* + itL]^* with the spin projection operator 
Pa (only for on-mass shell particles these two operators commute). In the present paper 
we take the spin-| propagator in the form [33] 

sr{p) = -^ ^p;^T* -pr(-p)> (2-13) 

where the spin projection operator is of the commonly adopted form in the Rarita- 

Schwinger formalism [40] 

1 2P^^P^ 1 
P^'{P) = g^^ - -YY - ^r^ - -^ il'-P" - YPn ■ (2.14) 

B. Fixing effective parameters 

The coupling constants, masses and cut-offs for the dipole form factors for the nucleon 
currents are taken the same as in the Bonn potential C [41], except for the u meson where 
the coupling is chosen as q^^nn = 11 (see also comments in Refs. [29, 30]). The choice 
of the magnitude of the coupling constant Qti'nn is still matter of debate. The OZI rule 
would imply a small value oi grf'NN-, while the above mentioned rj-rj' mixing conjecture can 
relate the corresponding coupling constants and express Qr/NN via Qr^NN providing in such 



a way a relatively large Qrf'NN [18, 30]. An analysis based on SU{3) symmetry and rj-r]' 
mixing angle, as performed in Ref. [30], has established an upper limit of Qrf'NN = 6.1. 
Note that in the analysis [30] some 50% of the pseudoscalar-pseudovector admixture in 
the NNrj' Lagrangian has been adopted. Subsequent investigations have shown that such 
a constant is too large, and a new upper limit Qrf'NN < 2 has been proposed [17] for the 
pseudo-vector coupling. Moreover, even values consistent with zero are also admitted 
in attempts to estimate the couplings and masses of possible excitations of the missing 
resonances [18, 31]. Evidently, the choice of a small value Qti'nn implies a neghgible 
contribution of the nucleonic currents, therefore increasing the role of meson conversion 
and nucleon-resonance currents in the rj' production. 

In our calculations we take Qtj'nn — 1-42 as reported in Ref. [2] and recently confirmed 
by the CLAS data [15]. The effective coupling constants for the resonance currents, 
whenever possible, have been estimated from the known decay widths of direct decay into 
1]' channels or radiative decay with subsequent use of the vector meson dominance (VMD) 
conjecture. The few remaining unknown cut-off parameters are taken either as constant 
or are adjusted to the experimental data. These parameters are listed in Tab. I. 

TABLE I: Resonance Parameters. For the spin-| resonance P13 the ofF-shell parameter is taken 
as z = —1/2 and the second coupling constant with vector mesons [29] is given in parenthesis. 





Su 


(1650) 


Pii 


(1710) 


^13 


(1720) 




qmnn* 


A [Gey] 


qmnn* 


A[GeV] 


Qmnn* 


k[GeV] 


IT 


1.47 


1.2 


1.47 


1.2 


0.2 


1.2 


V 


0.7 


1.2 


1.9 


1.2 


0.6 


1.2 


CO 


2.8 


1.2 


6.28 


1.2 


10(2) 


1.2 


p 


1.1 


1.2 


2.1 


1.2 


-10(-6) 


1.2 


v' 


1.18 


0.9 


1.4 


1.2 


1.0 


1.2 



All vertices with off-mass shell nucleons and nucleon resonances are dressed with (i.e. 



are to be multiplied by) a form factor 



A' 



The coupling constants Quju^ri' and gpp^^' for the meson conversion current have been 
derived from a combined analysis of radiative decays within the vector meson dominance 
(VMD) model and within effective Lagrangians with SU(3) symmetry providing 

^^.,'^4.90, (?,,,, ^5.84. (2.16) 

Note that naive direct calculations of these constants within the VMD model alone 
can provide slightly larger values, e.g. gvvri' ~ 6.5 — 7.0. The corresponding cut-offs, 
^Luujri' = ^pprf = 1-2 GeV , determine the form factor in the conversion vertex rj'VV (where 
V denotes u) or p) which is chosen as (see Refs. [27, 30]) 

Fyy,,{K, kl kl) = j,2_flj,2_k2 - (2-17) 

In accordance with the procedure of determining the coupling constant, the form factor 
(2.17) is normalized to unity if one vector meson is on-mass shell, while the other one 
becomes massless, e.g., Fvvrf{-^,k1 = rriy^kl = 0) = 1. In calculations performed in 
Ref. [31] the form factor (2.17) has been changed from the monopole to dipole form. 
As a result, the conversion contribution becomes essentially suppressed in comparison to 
previous results [30]. 

III. THE REACTION pp -^ pp7e+e 

For the Feynman diagrams in Fig. 1 generated by the Lagrangians (2.3)-(2.10), the 
invariant amplitude T can be cast in a factorized form 



The amplitude Tf^j^_^j^^ , describes the production process of a pseudoscalar meson (an 
off-mass shell r]' meson in the A^A^ collision), while the second amplitude T ,i^ +g- de- 
scribes the Dalitz decay of the produced meson into a real photon and a di-electron. In 
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the propagator of the ri meson the mass tti^/ is replaced by vn^i — iV^i /2 to take into 
account the finite hfe time of the r]' meson. 

For such factorized Feynman diagrams one can separate, in the cross section, the depen- 
dence on the variables connected with the Dalitz decay vertex and perform the integration 
analytically [26, 27]. Then, formally, Eq. (3.1) allows one to write the differential cross 
section (2.2) in a factorized form as well, 

da _ ai r)i^^e+e- 1 1 rl^^tot (o o\ 



where the integrals over the final di-electron and photon variables have been carried out 
analytically (see for details [26]). The decay rate dV /ds-f* for the rj' meson is defined as 

dVr^i^^e+e- _ '2^CXem /-, ^rf\ p .p . ^|2 /o oA 

The electromagnetic fine-structure constant is denoted as aem, and F^/^^. (s^.) is the 
transition rj' — > 77* form factor. The cross section for the production of a pseudoscalar 
meson with rj' quantum numbers but with s^' 7^ rn^, is 

j5 iot _ \ 

NN^NN^' 2(27r)yA(s,m2,m2) 

Xt E \T''NN^NN^'?dsi2dR2{NiN2 ^ ^^12)^^2(512 ^ N[N'^), (3.4) 

spins 

where S12 is the invariant mass of the two nucleons in the final state, and the two-particle 
invariant phase-space volume R2 reads 

fl,(„i ^ ed) = V^'^->-<."'S),i,.. ,3,5) 

8 Sab 

It can be seen from Eq. (3.2) that all the peculiarities of the cross section for the full 
reaction A^A^ -^ NN'~fe~^e~ are basically determined by the sub-reaction A^A^ -^ NNrj' 
with creation of a virtual rj' meson. Hence, before analyzing the full reaction, we shall 
proceed with a detailed study of the sub-reaction A^A^ — > NNt]', i.e. the production of 
an on-mass-shell r]' meson. 
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A. vj production in NN -^ NNr/' reactions: initial and final state interaction 

It ought to be mentioned that generally the Feynman diagrams deal with asymptoti- 
cally free particles, i.e. account for reactions with non-interacting particles in initial and 
final states. However, in the real process (2.1) the two nucleons can interact in the initial 
state (ISI) before the Vj creation, and in the final state (FSI) as well, provoking distortions 
of the initial and final A^A^ waves. 

The initial state distortion due to the A^A^ interaction before the rj creation is to 
be evaluated at relatively high energies, i.e., larger than the threshold of the Vj meson 
production, Tfej„ ~ 2.41 GeV . Therefore, one can expect that the variation with the 
kinetic energy of ISI effects is small. As shown in Ref. [42], the effect of ISI can be 
factorized in the total cross section and henceforth it plays a role of a reduction factor in 
each partial wave in the cross section. This reduction factor depends on the inelasticity 
and phase shifts of the partial waves at the considered energies. At the threshold, the 
number of initial partial wave is strongly limited by the partial waves of the final states 
and, in principle, one can restrict oneself to ^Pq and ^P\ waves. Experimentally it is 
found [43] that at kinetic energies of the order of few GeV the phase shifts ^Pq and ^Pi 
are indeed almost energy independent and the reduction factor for each partial wave can 
be taken constant. In our calculations we adopt, for the reduction factor (^, the expression 
from Ref. [42], which leads to (^ = 0.277 for the ^P\ wave and C^ = 0.243 for the '^Pq waves 
(cf. Ref. [29]). 

Final state interaction effects in the A^A^ system have been calculated within the Jost 
function formalism [44] which reproduces the singlet and triplet phase shifts at low ener- 
gies. Details of calculations of FSI with the Jost function can be found in Ref. [45]. 

B. Results for r]' meson production in NN collisions 

The amplitude T^j^^jy^^,, besides the above listed meson conversion contribution, 
includes the nucleonic and resonance currents, each of them being a coherent sum over 
all the considered exchanged mesons. The corresponding parameters are described above 
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in Section II B. It is worth mentioning that, in spite of the large number of considered 
diagrams and the large number of effective parameters, there is not too much freedom 
in fitting the cross section as long as the known resonances are considered. However, if 
one fits data with additional inclusion of missing resonances, a bulk of parameters remain 
practically free and can be fine-tuned to optimize the description of data. We performed 
some investigation of such a model with unknown resonances by considering few missing 
resonances with masses about 2GeV/(? and found that, since the energy dependence of 
the cross section is rather smooth, various sets of parameters fit equally well the data 
near the threshold. However, for each of the considered set it turns out that the cross 
section rises too rapidly at large values of the excess energy. Therefore, to reconcile this 
dependence with a smoother behavior at large excess energies one can, in principle, add 
more and more resonances and adjust appropriately the relative (unknown) signs of the 
couplings to compensate contributions from different diagrams (the interested reader can 
find a detailed analysis of such a situation in Ref. [31]). 

In our present calculations we included only the known four-stars resonances S'ii(1650), 
Pii(1710) and Pi3(1720). Results of numerical calculations for the total cross section 
^NN^NNri g^j,g presented in Figs. 2 and 3, for proton-proton and proton-neutron reac- 
tions, respectively. Experimental data are taken from [46]. The total contribution of 
nucleonic and resonance currents (dot-dashed lines) is found to be much smaller than the 
mesonic diagrams alone providing a constructive interference in pp and a destructive one 
in pn reactions. As already discussed, the small contribution of the nucleonic and reso- 
nance currents is due to the smallness of the coupling constants gr,'NN and gr,'N*N- The 
magnitude of other parameters entering in our calculations are dictated by the known 
experimental data. A reasonable agreement with data is achieved almost by the contribu- 
tion of the meson conversion currents. This does not contradict to the results reported in 
Ref. [30], however, if one adopts a dipole form factor [31] for the meson conversion vertex, 
then additional resonances with freely adjustable parameters have to be accommodated 
in the model. 

To further investigate the role of different diagrams, one needs to consider quantities 
being more sensible to the production mechanism. For instance, a detailed investigation of 
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the angular distribution (see Figs. 4 and 5) shows that the mesonic current provides a form 
of the differential cross section which evolves from a flat to a convex curve with increasing 
excess energy. Contrarily, the nucleonic and resonance currents give a concave shape of 
the corresponding contributions. The interference effects (which are maximum in forward- 
backward directions) lead to a slightly convex form of the resulting curves. Experimental 
data are still too scarce to determine more precisely the relative contributions of different 
diagrams. Note again that, taking into account more resonances, the description of the 
experimental data on angular distribution at larger excess energies cannot be improved 
[31]. New data will essentially enlighten the problem. 

Note also that, even achieving a good fit of the cross section in proton-proton reactions, 
it is not a priory clear whether the obtained set of parameters equally well describes also 
the proton-neutron reactions. The isospin dependence of the amplitude is determined by 
a subtle interplay of different diagrams with different exchange mesons (scalar, vector, 
iso-scalar, iso- vector etc.). Once the parameters for pp reactions are fixed, the pn ampli- 
tude follows directly from this set of parameters without any further readjustment. Of 
course, the ISI and FSI factors are different for pp and pn systems, but they are fixed by 
independent information. Our results displayed in Fig. 3-5 may serve as predictions for 
the pn channel. 

Figure 6 exhibits the isospin dependence of the total cross section. Similar to u;, and 
?7 production [26, 32] the ratio of the pn channel to the pp channel depends on the excess 
energy. This means that a simple (i.e. constant) isospin factor (see discussions in Ref. [47]) 
cannot relate these channels. In addition, our calculations point also to the possibility of 
rather different angular distributions at the same excess energy, as evidenced in Figs. 4 
and 5. The angular distribution in the pn channel is fairly fiat. The rise of the isospin 
dependence at low excess energies is entirely determined by the FSI effects which are 
rather different in the pp and pn channels. At larger excess energy the FSI effects vanish 
and the isospin ratio is determined entirely by the corresponding Feynman diagrams. 
Figure 6 exposes the importance of the pn channel for heavy-ion collisions. This seems 
to hold for the production of light vector mesons and virtual (e^e^) bremsstrahlung too 
[33]. 
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IV. DALTZ DECAY 



A. Preliminaries 



Consider first the general process of Dalitz decay of an on-mass shell pseudo-scalar 
meson P into a real and a virtual photon (di-electron). Within the present approach this 
reaction can be considered as a two-stage process, when the meson decays firstly into two 
photons, and secondly, one of the virtual photon decays into a di-electron pair. The decay 
width of production of two real photons is calculated from Eq. (2.10) as 



s 



3/2 
P r2 



fp,M (4.1) 



and serves for a determination of the coupling constant /p-y-y from experimental data. 
The square of the 77 invariant mass is denoted by sp. Since in what follows we are 
interested in production and decay of rj' mesons, being generally off-mass shell, we use 
the notation s^/ instead of Sp, bearing in mind that s^/ 7^ rn^,. Contrarily to the vector 
meson case [26], instead of the factor 1/3 (due to averaging over three projections of 
the spin of the vector particle) in Eq. (4.1) a factor of 1/2 appears due to two photons 
in the final state. Equation (4.1) yields \fr)''y'y\ — 0.0314 GeV^^ for the known width 
Fj,/-,^^ = (4.30 ± 0.15) keV [34]. In our calculations the sign of the coupling constant has 
been taken positively. 

The transition form factor -F^'^^* (S7*) requires separate consideration. As stressed 
above, the functional dependence of form factors upon the momentum transfer is a source 
of information on general characteristics of hadrons, such as charge and magnetic dis- 
tributions, size etc. Also, form factors are known as important quantities characterizing 
bound states within non-perturbative QCD. The main theoretical tools in studying ex- 
clusive processes within the non-perturbative QCD are approaches based on light cone 
sum rules and factorization theorem (see, e.g., [48-52] and references therein quoted). At 
high Q^, i.e., at high virtualities of quarks the hard gluon exchange is predicted to be 
dominant. However, since the virtuality of a quark is determined by Xip {xi is fraction 
part of the total momentum p carried by quark) the QCD asymptotic regime for low 
values of Xi can be postponed until extremely large values of the momentum transfer and 
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non-perturbative (soft) contributions play an important role in exclusive processes. In 
this connection an investigation of form factors in large intervals of momentum transfer, 
including the time-like region, serves as an important tool to provide additional infor- 
mation about the QCD regimes and the interplay between soft and hard contributions. 
Besides, there is a number of more phenomenological models, e.g., based on the dispersion 
relation technique [53, 54], or on use of vector meson dominance models (for details see 
Refs. [55-58]), or on effective SU(3) chiral Lagrangians with inclusion of the f/(l) non- 
Abelian anomaly [58-60] . It is also known that the asymptotic behavior of the transition 
form factors can be determined from the corresponding axial anomaly [61]. Traditionally, 
the electromagnetic form factors are studied by electron scattering from stable particles 
which provides information in the space-like region of momenta where, as well known, 
the experimental data can be peerless parameterized by dipole formulae. This in turn 
means that in the unphysical region, i.e., for kinematics outside the reach of experiments 
with on-mass shell particles, the analytically continued form factors display a pole struc- 
ture. Intensively studied form factors are the ones of the light pseudo-scalar mesons, in 
particular the pion. Heavier pseudo-scalar meson form factors have received less atten- 
tion since their experimental determination is more difficult. However, there are already 
experimental data on form factors of the heavier pseudoscalar mesons, such as rj and rj', 
most of them in the space-like region [25, 62]. With the spectrometer HADES [22] which 
can investigate with a high efficiency di-electrons production in pp collisions in a wide 
kinematical range of invariant masses, an additional study of the transition form factor 
for 1]' —^ 7e^e~ process becomes feasible. 

The simplest and quite successful theoretical description of form factors can be per- 
formed [55, 57, 58] with the VMD conjecture, within which a reasonably good description 
of elastic form factors in the time-like region has been accomplished. By using the current- 
field identity [55, 58] 

2 9 2 

777 7T7 777 / 

J' = -ej^Ko - ^t"^'^ - ^^K + ■ ■ ■ (4-2) 



with the coupling constants /-^p, f^^i and f^^ known [2, 55, 63, 64] from experimentally 
measured electromagnetic decay widths, one can also compute the transition form factor 
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FX^P{sy*) by evaluating the corresponding Feynman diagrams. The result is 

F;,^^F(v) = V4^^^ E /'''7 ."1 ^ E C^v^^^, (4.3) 

V=p,cj,(}!)... h'llJlV "^V -^7* V=p,cj,(}!)... ""^V -57* 

where the summation over vector mesons V , besides p, 00 and mesons, contains 
also heavier vector particles [55]. The normalization of the form factor at the origin, 
Fri'.y^*{sY* = 0) = 1 implies that J^Cy = 1. Actually, in real calculations one takes into 
account only a restricted number of vector mesons, for which J^Cy 7^ 1, i-e. in to obtain 

V 

correct normalization, the coefficients Cy must be rescaled. Here,we take into account 
p, uj and mesons. The absolute values of Cy in Eq. (4.3) have been calculated from the 
rates Eq. (4.1) as 

TiV ^ e-e-) = i^!^. (4.4) 

3 J^y 

The relative signs have to be fixed from some more general considerations. For instance, 
SU(3) symmetry implies an opposite sign for f^^ relative to the corresponding couplings 
for p and a;, viz. f^p : f^^ : /^^ = 1:3: (— 3/\/2). We suppose that the sign of 
the coupling constant f-q/y-y does not depend upon the mass of the vector particle, i.e., 
the ratio frfv-yl f-q'-y-y in the present paper is taken positive for all mesons. Numerical 
calculations of Cy with experimental decay rates [34] result in Cp = 0.946, C^^ = 0.095 
and C,j) = —0.041, with a rescaling factor ~ 0.83. The accordingly computed transition 
form factors (see also Refs. [5, 24, 65]) are in a reasonable agreement with the available 
experimental data [62], which usually are parameterized in a monopole form 

^.'77*(Q') = y3qya2;' (4-5) 

where the effective pole-mass parameter A^/ ~ O.SlGeV is a combination of the masses of 
the vector mesons. It is worth emphasizing that the experimental parameter A^/ has been 
obtained by fitting data mostly in the space-like region (Q^ < 0). In case of time-like 
arguments of the form factors, i.e., when Q^ = s^. > the pole mass parameters my 
in Eq. (4.3) receive a small imaginary part, proportional to the total V meson widths, 
my — ^ m,y—iry/2. For the monopole parameter A^/ in Eq. (4.5), we adopt the imaginary 
part as rp/3. 
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In Fig. 7 we compare the transition form factor computed by the VMD formula, 
Eq. (4.3), with the corresponding monopole experimental fit (4.5). It is seen that the 
form factors obtained within VMD model and fitted to data are in a good agreement in a 
wide interval of di-electron masses. Note that in the Dalitz decay Vj -^ ^e^e" the p and uj 
poles occur in the physical region allowed for the invariant mass of the di-electrons. This 
implies that an experimental investigation of the e+e^ mass spectrum in the vicinity of 
the p and uj poles will provide an excellent test of the validity of the VMD conjecture at 
large virtualities in the time-like region. 

V. THE COMPLETE NN -^ NN7e+e REACTION 

In order to emphasize the dependence upon the di-electron invariant mass we rewrite 
the cross section (3.2) in the form 

da 2aerrnr- I \>2 Br{r]' ^ 2-f) 



rr]'^^*{Sy*)\ 2 -^l-^^ 



c^V 37r2 1^"^^ ^-^^1 m^, ^ ^^' 



^(V) = / d^ ^^ /^' ^/" .y N^N^NNAO, (5.1) 



S,max 

^A^ 

where we introduced dimensionless variables a,b and ^ as follows 

W — ) (^ — , ^ — , qmax — W-^J 

and Br{ri' -^ 27) is the branching ratio of the rj' meson decay into two real photons, 
Br{ri' -^ 27) = (2.12 ± 0.14)% [34]. Since the total width of the rj' meson is fairly small, 
the parameter b in (5.1) provides a sharp maximum of the integrand function at ^ = a 
as far as the parameter a > 1 (or, equivalently, the di-electron invariant mass s* < mi,). 
This allows to take out from the integral the smooth function 0"/^7v->Ar7Vr/ ('C = '^) ^"^^ 
calculate it at s,,/ = ml,, i.e. at the rj' meson pole mass. The remaining integral can then 
be carried out analytically. However, at a < 1 when the di-electron mass is larger than 
the 7]' pole mass, i.e., the integrand does not exhibit anymore a resonant shape, and the 
integral needs to be calculated numerically. 

Equation (5.1) shows that the di-electron invariant mass distribution da/ds^* is propor- 
tional to the transition form factor F^/^^*(s^*) so that measurements of this distribution 
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provide a direct experimental information about this quantity. At values s^* < mf, the 



smooth part of the integrand, a '^ ** (0? can be pulled out the integral at 



X(V) = ' 2 ^^^^^^"(V = S')- (5-3) 



^ = a obtaining 

n(a-l? 

Equations (5.1) and (5.3) allow to define an experimentally measurable ratio which is 
directly proportional to the form factor 

^^^f'' (da/ds *)\ s* ■ \ l-s./m^, / IF. ,(s -* OP' ^^ 

where s^*min is the minimum value accessible experimentally (in the ideal case this is 
the kinematical limit SY*min = 4:fil). At low enough values of SY*min, the transition form 
factor is close to its normalization point F^'^^*(0) = 1 and the ratio (5.4) is just the 
transition form factor as a function of s^* . 

As a approaches unity keeping the maximum position still within the integration range, 
one can again take out the integral the smooth function q-^'^^'^^v ^^ g^, = .^2 ^ However, 
now the function {C, — 1)^/'^ cannot be considered smooth enough and must be kept under 
the integration. Nevertheless, even in this case the remaining integral can be computed 
analytically [27] and one can still define a ratio analogously to Eq. (5.4) which allows for 
an experimental investigation of the form factor near the free r]' threshold (s^* -^ 1^1') ■ 
At s^* > rn^i the integral X(s^*) does not have anymore a sharp maximum and it ought 
to be calculated numerically. 

1 /2 

In Figs. 8 and 9 the di-electron invariant mass distribution da/ds^* is exhibited as a 
function of the invariant mass Js^* calculated by (i) Eq. (5.1) with the transition form 
factor F,rf^^*{s^*) from the VMD model (solid line), (ii) with the monopole experimental fit 
(4.5) (dashed line), and (iii) without accounting for the form factor, i.e. with F^/^^. (s^*) = 
1, known also as QED calculations (dotted line). It is clearly seen that without form 
factors the calculated cross section rapidly vanishes as Js^* — > mj,/, while the inclusion 
of the form factors results in a resonant shape of the distribution. 

It should be noted that the considered final state (two nucleons, one real photon and one 
di-electron) can be produced by various competing channels too. For instance, prominent 
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channels are Dalitz decays of pions and 77 mesons as well. However, the contribution 
of different pseudo-scalar mesons to the di-electron mass distribution is well separated 
kinematically. Thus, at low values of Sy* the pion contribution is predominant [21]. At 
intermediate s^* , close to the 77 meson pole mass, the invariant contribution of rj Dalitz 
decay is much larger than the corresponding contribution from rj' mesons. This can be 
understood from an inspection of Eq. (5.1) which is valid also for rj meson production. The 
main difference in mass distributions occurs from different values of the total cross section 
(7*N^N->NNP fo^ V s-iid 7]' production. For an illustration, in Fig. 10 we present a comparison 
of the total r] and r]' cross sections as a function of the excess energy As-*^/^. It is seen that 
the shape of the energy dependence is basically the same, however, the 77 production cross 
section, at equal excess energy, is larger by more than one order of magnitude. Moreover, 
at equal initial kinetic beam energy the r] curve is essentially shifted towards higher excess 
energies. For instance, at T = 2.bGeV the excess energy for rj' is As^/^ ~ 32MeV^ while 
for 7] meson As^/^ ~ 442 MeV . In this case, the r] cross section becomes larger than the r]' 
cross section by more than three orders of magnitude. The remaining terms in Eq. (5.1) 
are of the same order for both r] and rj' mesons. Hence, at s^* close to the r] pole mass 
the invariant mass distribution, for the specified final state with one real photon and one 
di-electron, is entirely governed by the rj Dalitz decay. However, at larger di-electron 
masses, where the (virtual) r] meson is highly off-mass shell, this contribution rapidly 
vanishes [27] and, consequently, beyond the rj pole mass the distribution is governed by 
the 77' Dalitz decay. 

Nevertheless, some comments about this kinematical region are in order. Since the 77' 
mass is above the p and uo pole masses, the transition form factor and, correspondingly 
the mass distribution, display a sharp peak due to smallness of the u meson width. To 
get a feeling on the required detector resolution to identify experimentally such a peak, 
we simulate it by 

da 



dst' 



1/2 




(5.5) 

In Fig. 11 we present the corresponding invariant mass distributions computed by taking 
into account this schematic detector resolution. It is seen that already for a resolution 
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characterized by c? = 15 MeV the uj peak essentially decreases and practically disappears 
at lower resolutions, d ~ 30 MeV. 

VI. SUMMARY 

In summary we have analyzed the di-electron invariant-mass distribution from Dalitz 
decays of rj' mesons produced in pp and pn collisions at intermediate energies. The cor- 
responding cross section has been calculated within an effective meson-nucleon approach 
with parameters adjusted to describe the free vector and pseudo-scalar meson produc- 
tion [26, 32] in A^A^ reactions near the threshold. In particular, we found a suitable 
parametrization of the cross section for pp — ^ ppi]', which can be used for a prediction 
of the channel pn — > pnr]', which is important in heavy-ion collisions. We argue that by 
studying the invariant mass distribution of the final e"'"e^ system from rj' Dalitz decays in 
a large kinematical interval one can directly measure the transition form factor F^/^^*, e.g., 
in pp collisions. Such experiments are possible, for instance at HADES, CLAS, KEK-PS, 
and our results may serve as predictions for these forthcoming experiments. Experimental 
information on form factors is useful for testing predictions of hadronic quantities in the 
non-perturbative QCD domain. In particular, as mentioned in Introduction, the use of 
the quark flavor basis in a combined analysis of rj and rj' mesons with one mixing angle, 
but with topological effects connected with the U{1) axial anomaly kept, provides an 
analytical form of the transition form factors with two terms, one related to non-strange 
quarks another one being entirely expressed via contributions from strange components 
[24]. This is in full analogy with the vector meson dominance with the meson incor- 
porated. Also, a precise determination of the nucleon-nucleon-?]' coupling constant Qnnti' 
together with an investigation of the polarized deep inelastic lepton scattering, in connec- 
tion with the Goldberger-Treiman relation [10], can supply additional information on the 
flavor- singlet axial constant Ga{0) and on the gluon coupling QNNg, hence understanding 
of the origin of the " spin crisis" . 

We predict the transition form factor and the invariant mass distribution to exhibit 
maxima at di-electron masses close to the p and u poles, in contrast to a pure QED 
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calculation without strong interaction effects which provides a smooth behavior of the 
distribution in the whole kinematical range. 
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FIG. 1: Tree-level diagrams for the process NN -^ NN^e^e~ within an effective meson-nucleon 
theory, a) Dalitz Decay of 77'-mesons from nucleonic and resonance currents, b) Dalitz decay of 
r/' mesons from internal meson conversion. 
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FIG. 2: Total cross section for the pp — > ppr/' reaction as a function of the excess energy 
As^'^ = (y^ — Im — m^i). The dot-dashed line depicts the mesonic conversion contribution 
(MEG only), while the dashed curve represents the total contribution of nucleonic and resonances 
currents (Nucl.+Res.). The solid curve is the total contribution. Data are from Ref. [46]. 
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FIG. 3: The same as in Fig. 3 but for the pn — > pnr]' reaction. 
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FIG. 4: Angular distribution of ry' mesons in the center-of-mass system in pp — > pp-q' reactions 
at three different excess energies. Line codes as in Fig. 2. Data are from Ref. [46]. 
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FIG. 5: Angular distribution of r]' meson in the center-of-mass system in pn -^ purj' reactions. 
Line codes as in Fig. 2 
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FIG. 6: Ratio of the cross sections pn -^ pnr]' to pp -^ ppij' as a function of the excess energy. 
Line codes as in Fig. 2. 
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FIG. 7: Transition form factor ri — > 77* computed within the VMD model (sohd hne) in 
comparison with the monopole fit (dashed line) F^i^^* = 1/(1 — s*/A^,) with A^/ = (0.81 — iFp/S) 
[25, 62]. 
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FIG. 8: Invariant mass distribution of di-electrons from rj' Dalitz decay in the reaction pp = 
ppr]' -^ pp'ye'^e~ at proton kinetic energy Tp = 2.5 GeV (corresponding to the excess energy 
^^1/2 ^ 31.5 MeV). Solid line - calculations with the transition form factor within the vector 
meson dominance model, dashed line - transition form factor fit from Fig. 7, dotted line - pure 
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FIG. 9: Same as in Fig. 8, but at T^i^ = 3.5 GeV corresponding to As^/^ ~ 342 MeV . 
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FIG. 10: Comparison of total cross sections for rj and r]' production in pp collisions as a function 
of the excess energy. Details of calculations and references for data of the r/ meson cross section 
can be found in Ref. [27]. 
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FIG. 11: Invariant mass distribution weighted with a Gaussian corresponding to a finite mass 
resolution of the detector. Solid (dashed) line corresponds to the dispersion d = 30 MeV {d = 
IbMeV), while dotted line is without smearing. For T/jj^ = 3.5 GeV. 
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